MA018862 


A First-Order,  Worldwide,  Ionospheric, 
Time-Delay  Algorithm 

JOHN  A.  KIOBUCHAR 


Approved  lor  public  rolooso;  dfttri button  unlimltol 


IONOSPHERIC  PHYSICS  LABORATORY  PROJECT  4643 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 

HAHSCOM  AFB,  MASSACHUSETTS  01731 


AIR  FORCE  SYSTEMS  COMMAND,  USAF 


SECURITv  CLASSIFICATION  OF  THIS  «>AG'  n*t»  Entered) 


(C . 

V 


REPORT  DOCUMENTATION  PAGE 

WF/r  INSTRUCTIO.NS 
BF.FCni:  COMPLETING  FORM 

1.  REPORT  NUMBER  j 2.  GOVT  ACCESSION  NO. 

AFCRL-TR  -75-C302  j 

j PFCl^TKT"!  CATALOG  N*.'»»«f  P 

<K  TITLE  (mnd  Subtitle)  - 

-4^T\TE  Of  M.RGA7  * period  cove*co 
1 / 

/'XJ'IR ST -ORDER,  >VOR I. DIVIDE,  IONOSPHERIC 
/'TIME-DELAY  ALGORITHM 

6 PF.RiOPMiN.3  ORG.  REPORT  NUMBER 

AFSG,  No.  324  ^ 

L-AtfTTTDR(,«) 

)yj  John  A.  Klobuchar 

8 '"'->*■)  H ACT  OR  GRANT  NUMBER(«7 

9 PERFORMING  ORGANIZATION  NAME  AND  A009ESS 

Air  Force  Cambridge  Research  Laboratories 
Hanscom  AFB 
Massachusetts  01731 


U.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Cambridge  Research  Laboratory 
Hanscom  AFP 
Massachusetts 

U MONITORING  AGENCY  NAME  A AUUHfc??FftW* 


X 01731  V&AF-  VMS 


15a.  DECLASSIFICATION/ downgrading 

SCMEOULC 


16.  DISTRIBUTION  STATEMENT  (o  I ihl  • Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebmtrect  entered  in  Block  20.  It  different  irom  Report) 


f 

o ■ 

.f 


18.  SUPPLEMENTARY  NOTES 

TECH,  OTHER 


A 


;-:A 


;<\ 


19.  KEY  WOROS  (Continue  on  reveree  aide  it  necetamry  end  Identity  by  block  number) 

Total  Electron  Content 
Range  Error  Correction 
Ionospheric  Time  Delay 


\\  • > 

\ \S 


\y 


20  >f%S^TRACT  (Continue  on  revere*  etde  It  neceeaen-  end  Identity  by  block  number) 

A first-order  model  algorithm  designed  to  reduce  ionospheric  time-delay 
errors  by  approximately  50%  rins  on  a world-wide  basis  for  single  frequency 
users  of  the  Global  Positioning  System/NA VSTAR  is  described.  The  algorithm] 
was  designed  for  greatest  accuracy  during  times  of  day  when  ionospheric  time -I 
delay  errors  are  expected  to  be  largest.  The  algorithm  is  available  in  several 
options  from  one  which  requires  only  a monthly  mean  solar-flux  value  to  one 
which  requires  9 daily  or  monthly  update  coefficients.  Several  approximations 


DD  , X 7i  1473  coition  or  i nov  «s  is  obsolete 


Unclassified 


SECURITY  CL  ASSiFIC  ATtON  C F THIS  PAGE  (When  Dete  f'ntered) 


0// 


i 

s 

" 


i 

\ ! 


Unclassified 

SECURITY  CLASSIFICATION  of  this  PAC€fWi*»i  DM*  Ertttr**) 


Contents 

1.  INTRODUCTION  5 

2.  CHARACTERISTICS  OF  IONOSPHERIC  TIME  DELAY  6 

3.  DESIGN  OF  THE  ALGORITHM  6 

4.  MODEL  ALGORITHM  1 

5.  GEOMETRY  CALCULATION  REQUIREMENTS  9 

5. 1 Earth  Angle  9 

5.2  Ionospheric  Geodetic  Latitude  and  Longitude  11 

5.3  Conversion  from  Geodetic  to  Geomagnetic  Latitude  12 

5.4  Finding  Local  Time  12 

5.  5 The  Cbliquitj  Factor  14 

6.  ALGOPiTIJM  LATITUDE  COEFFICIENTS  15 

6.  1 The  DC  T»rm  15 

6.2  The  Amplitude  Term  18 

6.3  The  Phase  of  the  Diurnal  Period  18 

6.4  The  Diurnal  Period  18 

7.  OPTIONS  FOR  MODEL  INCORPORATION  INTO  THE  SYSTEM 

USER  COMPUTER  19 

7.  1 Option  1 19 

7.2  Option  2 19 

7.3  Option  3 . 19 

7.4  Option  4 20 

8.  DETERMINING  ALGORITHM  LATITUDE  COEFFICIENTS 

OPERATIONALLY  21 

9.  SUMMARY  OF  STEPS  IN  MODEL  ALGORITHM  21 


3 


Contents 


10.  MODEL  ACCURACY  23 

11.  CONCLUSIONS  AND  RECOMMENDATIONS  24 

Illustrations 

1.  Two  Examples  of  Ac-.ual  Monthly  Average  Time-Delay  Data  7 

2.  Approximation  to  Earth  Angle,  A 10 

3.  Approximations  to  Finding  Sub-Ionospheric  Coordinates  11 

4.  A Numeric  Representation  of  the  Error  of  the  Approximation 

to  Geomagnetic  Latitude  13 

5.  Approximation  to  Slant  Factor,  Z 14 

6.  Cosine  Model  Coefficients  (Bent  Model)  vs  Geomagnetic 

Latitude  for  Four  Seasons,  1969  16 

7.  Cosine  Model  Coefficients  (Bent  Model)  vs  Geomagnetic 

Latitude  for  the  Month  of  March  for  the  Solar  Minimum 
Year  1964  and  the  Solar  Maximum  Year  1969 

•* 


17 


A First-Order,  Worldwide, 
Ionospheric,  Time-Delay  Algorithm 


1.  INTRODUCTION 

The  development  of  an  advanced  satellite  navigation  system  having  worldwide 
coverage  and  positional  accuracy  of  a few  feet  has  been  undertaken  by  the  Defense 
Department  as  a Tri-Service  project  with  management  responsibilities  assigned 
to  the  Air  Force's  Space  and  Missile  Systems  Organization  (SAMSO).  This 
advanced  development,  called  the  Global  Positioning  System  (GPS)  or  NAVSTAR, 
will  use  timing  signals  a'  L band  emitted  from  satellites  in  12-hour,  60-degree 
inclination,  retrograde  orbits  to  determine  the  user  location.  Accuracies  of  a 
few  meters  rms  will  be  available  to  the  primary  class  of  users.  The  time  delay 
due  to  the  group  path  retardation  of  the  ea  th's  ionosphere  can  be  as  high  as  100 
meters  at  L band,  requiring  this  class  of  ui  irs  to  utilize  a two-frequency  receiv- 
ing system.  An  operational  two-frequency  system  with  two  L-band  signals  sepa- 
rated by  approximately  350  ts  400  MHz,  and  with  modest  receiver  signal -to -noise 
ratios,  is  expected  to  reduce  the  ionospheric  error  to  an  acceptable  value. 

There  will  be  another  class  of  users  of  the  GPS  who  will  not  require  the 
accuracies  attainable  by  the  two-frequency  system.  However,  they  will  still 
require  a first-order  correction  for  the  time  delay  of  the  earth's  ionosphere  which 
should  reduce  the  ionospheric  time  delay  error-  by  approximately  50%  or  more, 
on  an  rms  basis.  Due  to  the  user's  limited  computer  memory  and  computation 
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capability,  the  model  algorithm  must  be  of  minimum  size  and  complexity.  Such 
a first-order,  time-delay  correction  algorithm  is  the  subject  of  this  report. 


2.  CHARACTERISTICS  OF  IONOSPHERIC  TIME  DELAY 

The  ionospheric  time  delay  is  directly  proportional  to  the  total  electron  con- 
tent along  the  path  between  the  satellite  and  the  user.  Total  Electron  Content 
(TEC)  is  a function  of  many  variables;  among  them  are  local  time,  season,  geo- 
graphic location,  and  state  of  solar  and  magnetic  activity.  The  TEC  is  highest 
in  most  regions  of  the  world  within  a few  hours  of  local  noon  and  is  greatest  at 
regions  located  approximately  plus  and  minus  20  degrees  away  from  the  geomag- 
netic equator.  This  equatorial  anomaly  region,  as  it  is  called,  is  the  region  of 
the  world  where  single-frequency  GPS  users  will  encounter  the  greatest  errors. 
Statistics  of  daily  values  of  TEC  from  individual  stations  ranging  from  auroral  to 
equatorial  locations  show  that  the  rms  difference  from  the  monthly  mean  values 
during  the  midday  hours  is  approximately  20  to  30  percent.  Thus,  an  accurate 
specification  of  the  monthly  mean  values  alone  will  represent  an  elimination  of 
70  to  80  percent  (rms)  of  the  ionospheric  error.  The  problem,  then,  is  how  to 
best  represent  the  monthly  mean  behavior  of  ionospheric  time  delay,  particularly 
during  those  times  of  the  day  and  in  those  regions  of  the  world  when  the  time  delay 
is  greatest,  using  the  minimum  number  of  coefficients  and  computation  time. 


3.  DESIGN  OF  THE  LGORITHM 

The  model  algorithm  was  designed  to  meet  several  criteria.  First,  it  was 
designed  to  fit  best  during  those  times  of  the  day  when  the  ionospheric  delay  is  the 
largest,  that  is,  during  the  afternoon  period.  Second,  it  was  designed  to  have 
versatility  in  the  coefficient  array  so  that  trade  offs  can  be  made  in  the  number  of 
coefficients  used  to  represent  the  latitude  dependence  and  the  coefficients  can  be 
weighted  for  the  latitudes  of  greatest  interest,  such  as  the  Continental  United 
States  (CONUS).  Alternatively,  the  coefficients  can  be  weighted  for  best  fit  in 
those  geographic  locations  where  the  time  delays  are  greatest,  namely  the  regions 
located  approximately  plus  and  minus  20  degrees  either  side  of  the  geomagnetic 
equator.  Finally,  the  model  was  simplified  by  making  approximations  to  geometry, 
wherever  possible,  to  make  its  computer  running  time  shorter  and  storage-space 
requirements  smaller. 

The  available  experimental  values  of  TEC  from  which  time-delay  values  are 
obtained  are  not  sufficient,  by  themselves,  to  form  the  basis  of  a worldwide 


6 


time-delay  model.  However,  there  are  sufficient  TEC  data  to  yield  important 
information  on  the  statistical  behavior  of  time  delay  for  locations  representative 
of  several  regions  of  the  world.  These  statistics  will  be  presented  in  a separate 
report.  The  available  data  also  show,  quite  well,  the  diurnal  behavior  of  TEC 
for  various  seasons  and  solar-flux  values.  These  data  have  been  used  to  deter- 
mine the  form  of  the  d’urnal  algorithm  and  they  have  been  supplemented  with  the 
Bent1  model  representation  for  those  regions  where  actual  data  were  not  available, 
particularly  for  the  shape  of  the  latitude  gradient. 


4.  MODEL  ALGORITHM ' 


The  average  monthly  diurnal  behavior  of  time  delay  at  any  location,  as  a 
function  of  time  of  day,  has  been  represented  by  a simple  positive  cosine  wave 
dependence,  with  an  additional  constant  offset  which  we  shall  call  the  DC  term. 

In  Figure  1 two  examples  of  actual  ionospheric  time-delay  (T^)  data  are  shown 
with  a best-fit  cosine-type  curve  representation  of  each  set  of  actual  data  inclrded. 
While  the  cosine  curve  is,  of  course,  not  a perfect  fit  to  the  actual  data,  it  can  be 


Figure  1.  Two  Examples  of  Actual  Monthly  Average  Time-Delay  Data 


1.  Llewellyn,  S.  K. , and  Bent,  R.B.  (1973)  Documentation  and  Description  of  the 
Bent  Ionospheric  Model,  AFCRL-TR-73-0657  and  SAMSO-TR-73-252. 
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made  to  fit  reasonably  well  during  the  maximum  of  the  day  when  the  absolute  V 
values  are  greatest.  The  cosine  function  used  is: 


Tg  = DC  + A cos  . (1) 

where  DC,  A,  4,  and  P are  the  four  parameters  required  to  obtain  a complete 
diurnal  representation  of  Tg  at  any  location.  Only  the  positive  half  of  the  period 
of  the  cosine  curve  is  used  in  this  simple  representation;  thus,  when  the  absolute 
value  of  the  quantity  <t  - ^>)2v/P  becomes  greater  than  90  degrees,  Tg  is  repre- 
sented only  by  the  DC  term. 

The  simple  positive  cosine  representation  of  the  actual  Tg  is  designed  to  fit 

best  near  the  hours  of  maximum  T . The  fit  near  sunrise  and  well  after  sunset  is 

g 

expected  to  be  poor;  however,  during  these  times  the  absolute  value  of  Tg  is  much 
smaller  than  near  the  peak  of  the  diurnal  curve.  The  lower,  later  nighttime  values 
are  represented  only  by  the  constant  DC  term. 

A check  of  the  cosine  diurnal  fit,  as  compared  with  the  actual  monthly  mean 
values,  for  a few  months  representative  data  from  eight  different  TEC  observing 
stations,  revealed  that  the  cosine  algorithm  was  within  approximately  10  percent 
of  the  actual  24 -hr  rms  error. 

As  a further  computation  efficiency,  since  the  cosine  diurnal  dependence  is 
used  only  for  the  positive  half  of  the  waveform,  the  first  two  terms  of  the  cosine 
expansion  are  used  instead  of  the  actual  cosine  function  itself.  Thus,  the  algorithm 
becomes : 


T = DC  + A 
g 


L 2 24  J ' 


where 


_ (t  - 


(2) 


and  t is  local  time  at  the  ionospheric  point. 

The  error  between  the  expansion  and  the  actual  cosine  reaches  a maximum  of 
0. 02A  at  the  points  where  it  just  intersects  the  DC  value.  Since  the  cosine  fit  is 
too  low  at  these  points,  the  series  approximation  to  the  cosine  is  actually  slightly 
better  than  the  cosine. 

Other  forms  of  diurnal  algorithms,  such  as  a Gaussian  suggested  by  Antoniadis 
and  daRosa  (private  communication),  and  a statistical  F distribution  have  been 
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considered.  The  2-term  expansion  of  the  cosine  was  chosen  as  a reasonable  com- 
promise of  accuracy  versus  algorithm  difficulty. 


5.  GEOMETRY  CALCULATION  REQUIREMENTS 

In  determining  the  TEC,  ionospheric  time  delay,  from  a user  location  to  a 
satellite,  several  geometrical  calculations  must  be  made.  The  TEC  mu3t  be 
found  at  the  geographic  point  where  the  ray  path  intersects  the  mean  ionospheric 
height,  rather  than  at  the  user  location.  This  point  is  taken  here  at  a mean  height 
of  350  kilometers.  For  a satellite  at  the  lowest  elevation  considered  here,  5 
degrees,  this  point  is  approximately  14  degrees  of  earth  angle  away  from  the  user 
location.  The  obliquity,  or  slant  factor,  also  must  be  calculated  for  the  mean 
ionospheric  location.  Finally,  since  the  TEC  is  best  ordered  in  geomagnetic, 
rather  than  geodetic  coordinates,  the  assumption  is  made  that  the  TEC  is  only  a 
function  of  geomagnetic  latitude  and  local  time.  Any  two  points  on  the  earth  having 
the  same  geomagnetic  latitude  are  assumed  to  have  the  same  TEC  at  the  same 
local  time.  Thus,  a conversion  to  geomagnetic  latitude  is  required,  and  a calcu- 
lation is  made  of  local  time  at  the  geodetic  longitude  of  the  sub-ionospheric  loca- 
tion. If  the  exact  form  of  all  the  necessary  geometry  calculations  is  used,  the 
computer  time  used  would  be  greater  than  that  required  for  the  TEC  algorithm 
itself.  Therefore,  in  the  derivation  of  these  geometrical  functions,  simplifying 
assumptions  were  made  in  many  cases  to  reduce  the  calculation  complexity.  A 
discussion  of  each  of  the  steps  in  the  geometry  and  the  simplifying  assumptions 
follows. 

5.1  Earth  Angle 

The  ionospheric  time  delay  is  computed  at  the  point  where  the  ray  from  the 
satellite  intersects  the  ionosphere,  rather  than  at  the  observer.  For  an  observer 
located  near  the  earth's  surface  looking  at  a satellite  at  5 degrees  elevation,  the 
ionospheric  intersection  location  is  approximately  14  degrees  of  earth-centered 
angle  from  the  observer.  If  the  satellite  is  located  along  the  observer's  meridian 
at  5 degrees  elevation  angle,  the  14  degree  latitude  difference  can  represent  a 
significant  gradient,  particularly  near  the  equatorial  anomaly.  If  the  14  degree 
earth  angle  is  all  in  longitude,  at  40  degrees  latitude,  this  represents  a time  dif- 
ference of  1.2  hours,  during  which  time  the  ionospheric  time  delay  could  have 
changed  considerably. 

The  exact  earth  angle  representation  is: 


rw  : 


where 


Z = sin-1  (.948  cos  el)  , 

o 

for  a mean  ionospheric  height  of  350  km.  A useful  approximation  is: 


A s ,415-  _4 

A el  + 20  * 


In  Figure  2 both  the  actual  and  the  approximate  earth  angle  calculations  versus 
satellite  elevation  angle  are  given.  This  approximation  is  less  than  0. 2 degrees 
in  error  for  all  elevations  greater  than  10  degrees.  It  is  off  by  only  0.4  degrees 
and  0. 3 degrees  at  elevations  of  5 and  0 degrees,  respectively. 


A • 90-el-*in',[  cosell 


ELEVATION 


Figure  2.  Approximation  to  Earth  Angle,  A 


2.  R.  S.  Allen  assisted  in  developing  this  useful  approximate  form  of  earth  angle. 


5.2  Ionospheric  Geodetic  Latitude  and  Longitude 

Given  the  station  coordinates,  elevation  angle,  and  earth  angle  to  the  iono- 
spheric point,  its  coordinates  may  be  found  by  the  following: 


- -1 


<f>j  = sin  {sin  cos  A + cos  sin  A cos  az}  , 


XT  = \ + sin 
I o 


^-1 1 sin  A sin  az  f 
\ COS  (Jj  | 


(5) 

(6) 


where  is  the  station  latitude,  \Q  is  the  station  longitude.  If  a flat  earth  approx- 
imation is  used  the  following  simplification  can  be  made  to  the  above  two  equations: 


= <i>o  + A cos  az  , 


(7) 


xT  = x +AjsHL*L  . 

I o cos 


(8) 


In  Figure  3 the  geometry'  of  the  ionospheric  location  calculations  is  illustrated. 

At  latitudes  of  40  degrees  and  equatorward,  the  maximum  error  of  these  approxi- 
mations is  approximately  1 degree.  At  BO  degrees  latitude,  errors  in  determining 
ionospheric  latitude  reach  a maximum  of  -3  degrees  at  5 degrees  elevation  angle 


X,  = X0  + SIN'1  f SlN-A..SJ-N.-AZ-\ 
i COS  <f>z  1 


FLAT  EARTH  APPROXIMATION 


4>Z  " t 
XI  * K 


+ A COS  A2 
A SIN  AZ 
■ COS  <t> J 


Figure  3.  Approximations  to  Finding  Sub-Ionospheric  Coordinates 
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for  a few  azimuths.  Ionospheric  longitude  errors  reach  only  2 degrees  at  the 
same  location  and  direction.  -M  high  latitudes  above  75  degrees  an  ionospheric 
latitude  equal  to  the  observer  location  is  recommended,  regardless  of  observer'3 
elevation  angle.  Time  delays  at  these  latitudes  are  generally  much  lower  and 
more  highly  variable  than  in  the  equatorial  and  mid-latitude  regions.  Thus,  the 
time  delay  errors  introduced  by  neglecting  the  conversion  from  user  position  to 
mean  ionospheric  position  above  75  degrees  will  not  significantly  increase  the 
overall  worldwide  rms  errors  of  the  model. 

5.3  Conversion  from  Geodetic  to  Geomagnetic  Latitude 

Ionospheric  time  delay  is,  to  an  approximation,  best  ordered  by  geomagnetic, 
rather  than  geographic  latitude.  Thus,  it  is  necessary  to  perform  a transforma- 
tion from  geodetic  to  geomagnetic  latitude.  In  this  representation  there  is  no  true 
longitude  dependence  of  ionospheric  time  delay,  only  that  which  comes  about  due 
to  the  differences  between  geomagnetic  and  geodetic  latitudes  as  a functi-  n of 
longitude.  The  transformation  from  geodetic  to  geomagnetic  latitude,  assuming 

that  the  earth’s  magnetic  field  can  be  represented  by  an  earth  centered  dipole,  is 
v.  3 

given  by: 

sin  $ = sin  sin  $p  + cos  $ cos  $p  cos  (X  - Xp)  , (9) 

where 

*p  = 78.3-N  , 

Xp  = 291.  0°E  . 

The  approximation 

<&  = * + 11.6°  cos  (X  - 291)  (10) 

represents  the  exact  form  to  within  1 degree  at  all  geomagnetic  latitudes  equator- 
ward  of  ±40  degrees  geomagnetic  latitude.  At  geomagnetic  latitude  up  to  ±65 
degrees,  it  is  within  ±2  degrees.  Figure  4 is  a numeric  representation  of  the 
error  of  the  approximation  to  geomagnetic  latitude.  Over  the  entire  CONUS  region 
the  error  in  the  approximate  form  is  less  than  1 degree. 

5.4  Finding  Local  Time 

Given  the  universal  time  and  the  approximate  longitude  of  the  ionospheric 
point,  Xp  the  local  time  at  the  ionospheric  point  is  simply: 

3.  Davies,  K.  (1966)  Ionospheric  Radio  Propagation,  Dover  Publications,  New 
York.  ~~  " 


12 


I I « • I I ( 9 « I I I I • 


O' 

fu 

l I 

<02 
► r * 
^ o o 
r ^ x 

Vi 


•X  o 

con  cj 

UJO 

a <ssr 

« W l« 

«f  (a  ui  <r 

ou  — •-* 

w u of  •- 

3 w O W 

ZtOiu  < o 

M *0  L« 


iWMf.'H 
I I I t « 


4 « (UW  (U*<  H 
I I t I • I I 


J WWWMH^ 

• 9 • • 9 • t 


MMMM  *4«4 


w W^*< 
• • • • 


w-«  v-l  <H  4M 


«*  tv  i\i  fit  •*»  tr 


h « .';  N n r j 


ift 


««<MN  CMIO  1 e\ 

«c 


I I 


O • *r 
►* 

o o 
HO*'* 
Of  0. 

O U. 
; -i 

u*  * 

o « 

OKUJ 
UK  3 
O -I  -» 


lf\ 


I I I I I I 


h 9 *■  * r* 


»-i  M rv»  fo  if' 


> MMfl  M «4  *4 


h m ry  n 4 o 


13 


Figure  4.  A Numeric  Representation  of  the  Error  of  the  Approximation  to  Geomagnetic  Latitude 


(hours) 


(11) 


If  t is  greater  than  24,  make  t = t - 24  hours. 

5.5  The  Obliquity  Factor 

The  vertical  time  delay  at  the  sub-ionospheric  point  must  be  multiplied  by  an 
obliquity  factory  defined  as  the  secant  of  the  zenith  angle  at  the  mean  ionospheric 
height.  An  average  ionospheric  height  of  350  km  is  assumed.  The  exact  obliquity 
or  slant  factor  is: 

SF  = sec  {sin-1[.  948  os  el]  } . (12) 


An  approximation  which  is  within  2 percent  of  the  exact  value  for  elevation  angles 
from  5°  to  90°  is: 

SF  = 1 + 2 . (13) 

The  approximation  saves  o trigonometric  functions  at  the  expense  of  finding  the 
cube  of  a number.  Both  the  exact  and  the  approximate  slant  factors  are  shown  in 
Figure  5. 


14 


"."w  •,  'iiv  «j*  'j*.'w-i<L»Ti 


6.  ALGORITHM  LATITUDE  COEFFICIENTS 

The  four  coefficients  used  in  the  time  delay  algorithm  are  all  functions  of 
latitude,  season,  solar  flux,  and  magnetic  activity.  A complete  representation  of 
all  these  variables  in  each  of  the  coefficients  would  be  complex  indeed,  if  even 
possible.  What  has  been  done  here  is  to  address  only  the  monthly  average  latitude 
dependence  of  each  of  these  four  coefficients  for  various  seasons  of  a solar  max- 
imum year  when  the  highest  ionospheric  time  delays  are  expected.  Possible 
methods  of  determining  the  daily  variation  of  the  algorithm  coefficients,  or  of 
simply  including  a daily  update  factor,  are  discussed  in  a later  section  of  this 
report. 

6.1  The  DC  Terra 

The  DC  term  in  the  model  algorithm  represents  the  nighttime  behavior  of  the 
ionosphere.  Recent  evidence  from  combined  Faraday  and  group-delay  measure- 
ments from  signals  transmitted  from  the  geostationary  satellite,  ATS-6,  shows 
that  the  Faraday  data,  which  is  really  representative  only  of  TEC  up  to  approxi- 
mately 2000  km  of  vertical  height,  is  considerably  in  error  in  the  nighttime.  The 
Bent  model  also  was  integrated  in  vertical  height  onlj  up  to  2000  km.  The  new 
group-delay  data  available  thus  far  show  a nearly  constant  additional  TEC  above 
the  TEC  normally  measured  by  the  Faraday  effect.  Thus,  with  this  model  repre- 
sentation it  is  a simple  matter  to  increase  the  DC  term  to  account  for  the  additional 
electrons  not  measurable  by  the  Faraday  effect. 

The  magnitude  of  the  DC  term  versus  latitude  is  shown  in  Figures  6 and  7. 

Since  the  DC  term  is  generally  much  smaller  than  the  amplitude  term.  A,  it  may 
be  sufficient  to  represent  its  latitude  dependence  by  a simple  linear  dependence 
with  geomagnetic  latitude  either  side  of  a reference  latitude,  or  by  a single  con- 
stant for  all  latitudes.  Alternatively,  it  can  be  represented  by  a second-order 
polynomial  in  geomagnetic  latitude  covering  the  entire  latitude  range.  The  largest 
errors  will  occur  at  the  equatorial  anomaly  locations  with  a second-degree  poly- 
nomial fit  because  of  the  inability  of  a second-degree  equation  to  fit  the  complex 
shape  of  the  actual  DC  term  versus  latitude. 

In  the  actual  finding  of  the  coefficients  for  the  latitude  dependence  of  the  DC 
and  the  amplitude  terms  it  is  important  to  make  certain  that  they  do  not  yield 
negative  values  at  certain  latitudes.  A least  squares  polynomial  fit  to  a set  of 
experimental  points  can  yield  unexpected  values  in  portions  of  the  data  set.  The 
responsibility  for  ensuring  that  the  latitude  coefficients  remain  positive  can  be 
done  well  in  advance,  so  that  thoroughly  tested  coefficients  will  be  loaded  into  the 
user  computer. 
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6.2  The  \uplitude  Tens 

The  amplitude  of  the  cosine  fit  is  the  most  important  term  in  the  algorithm, 
and  is  the  most  highly  variable  with  latitude,  season,  and  solar  flux.  It  is  the 
parameter  that  should  be  changed  on  a daily  update  basis,  if  this  is  deemed  neces- 
sary and  practical  with  the  operational  satellite  system.  To  closely  represent  the 
amplitude  term  as  a function  of  geomagnetic  latitude  several  coefficients  are 
required.  From  Figures  6 and  7 the  amplitude  dependence  with  geomagnetic  lati- 
tude can  be  seen.  Note  that  the  CONUS  region  extends  approximately  from  35°N 
to  55°N  geomagnetic  latitude.  It  is  obvious  from  the  curves  in  Figures  6 and  7 that 
no  single,  fixed  shape  will  match  the  latitude  dependence  of  the  amplitude  term. 
Therefore,  it  is  necessary  to  use  some  form  of  polynomial  representation  for  the 
latitude  de  tendence  of  amplitude.  Because  of  the  numerous  inflection  points  in  the 
curve  for  some  seasons  a polynomial  of  degree  5 or  6 is  necessary  to  obtain  a 
close  fit.  However,  if  the  equatorial  minimum  in  the  amplitude  term,  which  is  par- 
ticularly obvious  during  the  equinoxes,  is  smoothed  over,  a second -degree  poly- 
nomial can  do  a fair  job  of  representing  the  amplitude  behavior  versus  geomagnetic 
latitude.  The  same  caution  regarding  negative  values  at  high  latitudes  must  be 
observed  for  the  testing  of  the  amplitude  coefficients  as  was  given  for  the  DC  term 
latitude  coefficients. 

6.3  The  Phase  of  the  Diurnal  \ ariation 

Happily,  the  phase  term  in  the  ionospheric  lime-delay  algorithm  is  not  a highly 
variable  function.  A plot  of  the  phase  term  for  a northern  midlatitude  station,  for 
all  the  monthly  averages  of  a six-year  period  covering  from  solar  maximum  to 
well  down  towards  the  solar  minimum,  revealed  that  the  phase  of  the  diurnal  max- 
imum rarely  changed  more  than  plus  and  minus  one  hour  from  a mean  of  14  hours 
local  time.  Results  from  the  Bent  model,  in  Figures  6 and  7,  and  from  actual 
monthly  mean  TEC  data  from  several  stations,  show  that  the  phase  lies  between 
12.  5 and  15.  5 for  most  of  the  latitude  range  covering  the  mid  and  equatorial 
latitudes.  Therefore  a constant  value  of  14  hours  local  time  was  taken  for  the 
phase  of  the  diurnal  variation.  If  a particular  sub-model  is  required  for  a limited 
geographic  region,  perhaps  some  improvement  can  be  made  by  making  the  phase 
a function  of  time  of  the  year;  but,  in  this  attempt  at  a worldwide  representation, 
a constant  value  for  all  latitudes  and  seasons  is  deemed  adequate. 

6.4  The  Diurnal  Period 

The  period  of  the  diurnal  cosine  representation  is,  of  course,  not  fixed  at 
24  hours.  In  fact  it  is  nearly  always  significantly  higher  than  24  hours  and  is  a 
function  of  geomagnetic  latitude  and  season  as  shown  in  Figures  6 and  7.  As  is  the 
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case  with  the  amplitude  and  DC  terms,  a second-degree  polynomial  with  geomag- 
netic latitude  is  recommended  as  the  function  to  represent  the  latitude  dependence 
of  the  cosine  period  term.  Again,  the  equatorial  anomaly  region  will  not  be  per- 
fectly represented  by  a second  degree  polynomial;  but,  the  fit  should  be  adequate. 


7.  OPTIONS  FOB  MODEL  INCORPORVTION  INTO  THE  SYSTEM  L'SER  COMPUTER 

There  are  several  optional  methods  by  which  the  operational  user  computer 
and  satellite  data-link  combination  can  determine  an  ionospheric  correction.  Four 
options  are  suggested  here. 

7.1  Option  I 

The  most  sophisticated,  and  also  the  most  costly  in  terms  of  coefficients  to 
be  transmitted  down  the  satellite  to  user  link,  is  to  transmit,  and  update  once 
every  24  hours,  3,  second-degree  polynomial  coefficients  for  each  of  the  three 
algorithm  coefficients,  DC,  amplitude,  and  period,  making  a total  of  9 coefficients 
transmitted  and  updated  each  day.  Eight -tit  accuracy  for  each  coefficient  is  more 
than  adequate,  making  a total  of  72  bits  of  ionospheric  information  in  each  message. 

7.2  Option  2 

One  step  removed  in  complexity  is  to  update  only  the  three  amplitude  coeffi- 
cients on  a daily  basis  and  to  have  the  user  computer  store  all  the  coefficients  for 
the  seasonal  and  solar-flux  dependence  of  the  DC  and  the  period  terms.  If  these 
coefficients  are  changed  monthly  and  are  given  for  two  values  of  solar  flux  the  user 
computer  must  store  12  times  2 times  6 coefficients  or  144,  eight-bit  numbers  on 
a permanent  basis.  The  user  computer  would  interpolate  between  the  two  sets  of 
stored  coefficients  depending  on  the  actual  solar  flux.  In  addition,  three  amplitude- 
coefficient  numbers  plus  a solar  flux  must  be  transmitted  down  the  satellite  to  user 
link.  The  link  will  carry  only  4,  eight -bit  numbers,  or  32  bits  of  ionospheric 
information  in  each  message.  Using  this  option,  the  user  still  has  the  capability 
of  obtaining  updated  amplitude  coefficients  and  the  master  station  computer  has  the 
option  of  weighting  the  daily  values  of  amplitude  coefficients  for  best  accuracy  over 
any  particular  region,  such  as  the  CONUS. 

7.3  Option  3 

A third  option  would  have  all  the  coefficients  permanently  stored  in  the  user 
computer  and  would  have  only  a daily  value  of  solar  activity  transmitted  down  the 
satellite  to  user  link.  With  this  option,  the  data-line  requirements  are  minimal. 
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but  the  user  computer  must  permanently  store  12  times  9 times  2,  or  216,  eight- 
bit  numbers.  With  the  current  improvements  in  Read  Only  Memories,  ROM's,  a 
256  X 8 chip  could  conveniently  hold  all  the  necessary  coefficients. 

7.4  Option  4 

An  extremely  simple  representation  of  the  DC,  amplitude,  and  period  functions 
of  geomagnetic  latitude,  whmh  would  not  require  an  ROM,  could  be  used.  In  the 
case  of  DC,  this  would  be  a constant  which  would  vary  only  with  solar  flux,  not 
with  season.  For  amplitude,  only  linear  gradients  from  the  geomagnetic  equator 
would  be  used.  For  the  period  term,  a constant  for  all  seasons  and  solar-flux 
values  would  be  used.  The  update  information  required  would  be  only  solar  flux. 
Thus,  only  one  number  would  be  specified,  that  is,  the  solar  flux,  or  some  similar 
number  which  would  be  a multiplicative  factor  on  the  amplitude.  With  this  option, 
the  shape  of  the  latitude  gradient  would  remain  fixed,  or  would  change  in  a simple 
manner  with  seasonal  solar  flux,  and  only  one  dimensional  variations  would  be 
allowed,  as  only  one  parameter,  namely  the  solar  flux,  is  specified. 

With  options  1 and  2 above,  a network  of  ionospheric  monitoring  stations  could 
be  used  to  determine  the  latitude  coefficients  of  the  algorithm  coefficients.  With 
options  3 and  4,  only  a daily  solar-flux  value,  easily  obtained  from  the  ii.r 
Weather  Service's  Global  Weather  Central,  would  be  required. 

There  are  other  options,  for  example,  transmitting  two  sets  of  amplitude 
coefficients,  one  for  users  in  the  CONXJS  region,  and  one  for  users  in  the  equatorial 
region,  or  separate  coefficients  for  each  hemisphere,  but  these  have  not  been  con- 
sidered in  any  detail. 


Table  1.  Summary  of  Options  Considered 


Variables  Transmitted 

User  Storage 
Requirements 

Bits  Transmitted  Per 
Message  '8  Bits  Per 
Number  Assumed) 

Option  1 

DC  (3).  Amp.  (3), 
Period  (3).  Total  = 9 

Minimal 

72 

Option  2 

Amp.  (3),  Solar 
Flux.  Total  = 4 

144  (8 -bit 
numbers) 

32 

Option  3 

Solar  Flux  Only 
Total  = 1 

216  (8 -bit 
numbers) 

8 

Option  4 

Solar  Flux  Only 
Total  = 1 

Minimal 

8 
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8.  DETERMINING  ALGORITHM  LATITUDE  COEFFiaENTS  OPERATIONALLY 

If  consideration  is  given  to  an  option  which  would  require  ionospheric  meas- 
urements to  determine  updated  coefficients,  this  can  be  done  with  as  few  as  six 
monitoring  stations  located  alor  “ the  American  eastern-longitude  meridian,  A 
single  station,  looking  at  satellites  at  5 degrees  of  elevation  angle,  can  make 
measurements  of  the  equivalent  vertical  time  delay  at  latitudes  up  to  plus  and 
minus  14  degrees  from  the  station  latitude.  Thus,  3 stations,  spaced  at  approxi- 
mately 30  degree  latitude  intervals,  can  be  used  to  monitor  the  latitude  gradients 
of  time  delay  over  the  entire  Northern  Hemisphere  along  one  longitude  sector.  A 
total  of  6 stations  could  monitor  the  entire  latitude  gradient  for  both  hemispheres 
along  one  meridian. 

Since  no  true  longitude  dependence  is  included  in  the  model,  measurements 
along  one  meridian  are  sufficient  to  specify  the  entire  global  representation. 

Stations  in  the  northern  hemisphere  might  include  one  under  the  equatorial  anomaly 
at  the  Panama  Canal  Zone,  one  in  the  mid-latitudes  near  Washington,  D.  C. , and 
one  near  the  auroral  latitudes  located  at  either  Goose  Bay  or  Great  Whale  River, 
Canada.  In  the  Southern  hemisphere,  stations  might  be  located  at  the  geomagnetic 
equator  at  Huancayo,  Pem,  near  the  Southern  equatorial  anomaly  peak  at  the 
latitude  of  Santiago,  Chile  and  perhaps  in  the  Falkland  Islands. 

The  data  taken  at  these  monitoring  stations  could  be  used  to  make  daily  values 
of  coefficients  of  latitude  gradients  of  the  parameters  which  go  into  the  algorithm, 
namely  A,  DC,  phase,  and  period.  The  most  important  coefficient  to  attempt  to 
determine  on  a daily  basis  is  the  amplitude,  A.  The  amplitude  coefficient  is  the 
most  important  and  also  probably  the  most  highly  variable  one  on  a day  to  day 
basis,  during  most  seasons.  If  any  daily  updated  latitude  coefficients  are  to  be 
computed,  tt*e  amplitude  representation  should  be  considered  first.  If  only  a single, 
general  factor,  perhaps  dependent  upon  the  solar  ionizing  flux  or  state  of  geomag- 
netic activity,  is  to  be  transmitted,  it  can  be  used  to  alter  the  amplitudes  at  all 
latitudes  equally  on  a percentage  basis.  In  this  case  no  ionospheric  monitoring 
stations  will  be  required.  Provision  is  made  for  this  factor  in  the  algorithm. 


9.  SLMMARA  OF  STEPS  IN  MODEL  ALGORITHM 


Given:  (bQ,  X.  , az,  el,  UT,  Julian  date) 
Find  earth  angle: 


A = 


445 

el  + 20 


- 4 . 
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Find  sub -ionospheric  latitude: 


^ = tj>o  + A cos  (az)  . 

If  is  greater  than  75°,  let  = 75°;  if  is  less  than  -75°,  let  = -75° 
Find  sub-ionospheric  longitude: 


\r  = X.  + 


A sin  (az) 


I o cos 


’I 


Find  geomagnetic  latitude: 


+ 11.  6 cos  (\j  - 291)  . 


Find  local  time: 

1 = Ts  + UT  • 

Find  slant  factor: 


SF  = 1 + 2 


Also  given;  solar  flux,  coefficients  of  latitude  dependence  of  algorithm  variables. 

Calculate  DC  - Dq  + Dj  + D2  <t>*  + D3  <3>j  . 

2 3 

Ampl.  = Aq  + A j $ j + A2  + Ag  . 

Period  = PQ  + P2  *j+  P3*j  . 


Note;  In  option  4 these  polynomials  are  replaced  by  the  following: 


DC  = 2 


. 3*  filial 

90 


Ampl.  = 19  -1*1  + ~ ~go— ~ 13  + 5 cos  ( 


m = month  of  the  year 
Period  = 32  hours 
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Next,  calculate: 


x = ( Period  ) 2ff  * 

Test:  If  absolute  x is  greater  than  1.  57, 

T = SF  * DC  . 
g 

End. 

If  absolute  x is  less  than  1.  57, 

Tg  = SF*  j DC  + Ampl.  [l  \ ‘ 

End. 


10.  MODEL  ACCURACY 

This  is  the  toughest  part.  How  much  increased  accuracy  do  you  buy  for  adding 
additional  coefficients  ? Particularly,  how  much  greater  accuracy  is  obtained  by 
actually  setting  up  a six -station  ionospheric  monitoring  network  and  sending  down 
various  numbers  of  coefficients  to  the  users  via  the  satellite  to  user  data  link? 

Also,  since  these  coefficients  are  probably  24  hours  old,  how  much  improvement 
do  you  attain  over  monthly  mean  corrections?  Finally,  even  if  it  is  decided  to  send 
just  one  number  down  the  satellite  to  user  data  link,  and  store  all  the  other  coeffi- 
cients in  a user  computer  ROM,  how  big  does  the  ROM  have  to  be  to  attain  what 
correction?  Is  there  a "knee  in  the  curve"  of  number  of  coefficients  versus  accu- 
racy’ 

Unfortunately,  the  answers  to  the  above  questions  are  largely  unknown. 

Limited  tests  of  the  algorithm,  using  second-degree  polynomials  for  the  geomag- 
netic dependence  for  the  DC,  amplitude,  and  the  period  terms,  have  been  made. 
These  values  were  compared  with  the  monthly  average  actual  behavior  of  TEC  at 
a few,  selection  stations  and  a better  than  50  percent  rms  average  correction 
capability  was  achieved.  More  testing  remains  to  be  done. 

The  problem  in  determining  the  polynomial  coefficients  is  that  the  best  present 
set  of  them  is  probably  available  from  the  Bent  Ionospheric  Model:  but,  little  actual 
TEC  data  are  currently  available  t,hat  show  the  accuracy  of  the  Bent  model  in  the 
important  equatorial  anomaly  regions.  As  the  GPS  becomes  operational,  better 
coefficients  could  be  found  from  a network  of  ionospheric  monitoring  stations  and 
these  new  coefficients  could  be  put,  in  retrofit  fashion,  into  each  user  computer  to 
replace  the  original  ones  that  can  be  determined  from  the  Bent  model. 
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11.  CONCLUSIONS  AND  RECOMMENDATIONS 


Various  trade  offs  in  number  of  model  coefficients  have  been  suggested. 

The  diurnal  behavior  of  ionospheric  time  delay  at  any  location  can  be  repre- 
sented adequately  by  a positive  half  of  a cosine  wave  with  a DC  constant  term  for 
the  nighttime  behavior.  The  amplitude  and  period  of  the  cosine  term  are  functions 
of  geomagnetic  latitude,  as  is  the  DC  term.  Since  the  geomagnetic  dependence  of 
these  parameters  changes  with  season,  they  are  represented  as  polynomial  func- 
tions with  different  coefficients  for  each  season.  The  simplest  form  of  latitude 
representation  of  the  algorithm  parameters  is  a constant  for  the  DC  term,  chang- 
ing only  with  solar  flux,  with  constants  for  the  period  and  phase  terms.  The 
amplitude  term  is  the  big  unknown  at  this  time  due  to  the  unknown  shape  of  the 
latitude  gradient  and  its  changes  with  season  and  solar  flux.  Therefore,  some 
form  of  model  representation,  which  allows  the  latitudinal  shape  to  be  changed  as 
a function  of  time,  at  least  seasonally,  is  recommended.  A second-degree  poly- 
nomial of  amplitude,  with  geomagnetic  latitude,  is  recommended  as  the  coefficients 
can  be  weighted  for  a particular  latitude  range  of  interest,  such  as  the  CONUS. 

More  testing  of  the  suggested  models  against  the  monthly  average  TEC  data 
base  should  yield  a better  idea  of  the  rms  improvement  versus  model  complexity. 
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